Probability densities and preservation of randomness in wave turbulence 
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Abstract 

Turbulence closure for the weakly nonlinear stochastic waves requires, besides weak nonlinearity, randomness in 
both the phases and the amplitudes of the Fourier modes. This randomness, once present initially, must remain over 
the nonlinear evolution time. Finding out to what extent is this true is the main goal of the present Letter. For this 
analysis we derive an evolution equation for the full probability density function (PDF) of the wave field. We will 
show that, for any statistics of the amplitudes, phases tend to stay random if they were random initially. If in addition 
the initial amplitudes are independent variables they will remain independent in a coarse-grained sense, i.e. when 
considered in small subsets which are much less than the total set of modes. 

Introduction 

The theory of stochastic wavefields in weakly nonlinear dispersive media has a long and exciting history which started in 
1 1929 when Peierls derived his kinetic equation for phonons in solids [Q. Applications of these ideas appeared in the physics 
F^ f the ocean and atmosphere |21 El EI El El El El j laboratory and astrophysical plasmas ^imEli Bose condensates and 
C^jionlinear optics , anharmonic crystals ^ El E| ■ Any attempt to give a fair historical review would be doomed in 
^lich a short letter and we refer an interested reader for further references to the book ^7] and a more recent review . 
• The common name that has arisen for all these approaches is Wave Turbulence (WT). 
. ^ ' WT closure requires, besides weak nonlinearity, randomness in both the phases and the amplitudes of the Fourier 
^^nodes. Namely, all the phases and all the amplitudes must be statistically independent of each other, in some sense, 
5-arid the phases must be uniformly distributed. Such an approach was recently formulated in 19, 21 as a generalization 
. ^f the Random Phase Approximation (RPA) much loved by the physicists which, in its traditional form, ignores the 
amplitude randomness ^7] . We even kept the same acronym RPA but now read it as "Random Phases and Amplitudes" . 
Below, in section l?m we define explicitly what we mean by RPA. RPA does not fix the shape of the probability densities 
of the individual mode amplitudes and, therefore, it allows one to consider wavefields with non-decaying correlations 
which is helpful because such long correlations tend to arise naturally in WT systems. In we used RPA to describe 
the arbitrary-order moments of the wave amplitude, and in |21| we extended this approach to describing the one-mode 
probability density function (PDF) and considered solutions for this PDF corresponding to intermittency. In these works, 
however, RPA was assumed (but not proven) to hold over the nonlinear time. 

Such a proof is the main goal of the present paper. We shall consider initial fields of the RPA type, and we will prove 
that the RPA properties are preserved (i.e. no phase or amplitude correlations are generated with accuracy sufficient 
for the WT closure) over the nonlinear evolution time. In order to do this we shall derive an evolution equation for the 
full multi-mode PDF which will turn out to be the Zaslavski-Sagdeev (ZS) equation J3] (a WT cousin of the Brout- 
Prigogine equation for anharmonic crystals |15[I16| ). We will show that, for any statistics of the amplitudes, phases tend 
to stay random if they were so initially. If, in addition, the initial amplitudes are independent variables they will remain 
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independent in a coarse-grained sense, i.e. when considered in small subsets which are much smaller than the total set of 
modes. 

The original paper by ZS ^31 was also devoted to the study of the applicability of the WT closure and, therefore, it is 
appropriate here to mention in which way our approach is different. First, ZS consider the nonlinear interaction arising 
from the potential energy only (i.e. the interaction Hamiltonian involves coordinates but not momenta). This restriction 
leaves out the capillary water waves, Alfven, internal and Rossby waves, as well as many other interesting WT systems. In 
our work we remove this restriction by considering the most general three-wave Hamiltonian equation and we show 
that the multi-mode PDF still obeys the ZS equation in this case. Secondly, ZS studied the phase statistics only, whereas 
our work considers both the phases and the amplitudes because the amplitude statistics is as important for the RPA 
closure as the phase statistics. Thirdly, ZS presented an argument that the nonlinear frequency correction removes the 
need for the initial phase randomness, whereas we only state the preservation of the initial phase randomness. However, 
the ZS criterion for phase randomization was obtained from a rather non-rigorous (although highly intuitive) physical 
argument whereas our results follow from a systematic asymptotic expansion outlined in this Letter and the details of 
which will be published in a more extended paper j20) . 

The validation of the RPA properties gives this technique the status of a well-justified approach which, due to the 
simplicity of its premises, is a winning tool for the future theory of non-Gaussianity of WT, its intermittency and 
interactions with coherent structures. 

2 Statistical setup. 

Let us consider a wavefield a(x, t) in a periodic cube of with side L and let the Fourier transform of this field be ai{t) 
where index l^Z^^ marks the mode with wavenumber ki = 2Trl/L on the grid in the d-dimensional Fourier space. For 
simplicity let us assume that there is a maximum wavenumber k^ax (fixed e.g. by dissipation) so that no modes with 
wavenumbers greater than this maximum value can be excited. In this case, the total number of modes is TV = {kmax/'^L)'^ . 
Correspondingly, index I will only take values in a finite box, I G Bn C Z'^ which is centered at and all sides of which 
are equal to kmax/i^L — N^^^. To consider homogeneous turbulence, the large box limit iV — > oo will have to be taken. ^ 
Let us write the complex a/ as a; = Ai^i where Ai is a real positive amplitude and ipi is a phase factor which takes 
values on 5^, a unit circle centered at zero in the complex plane. Let us define the iV-particle joint PDF V^^^ as the 
probability for the wave intensities Af to be in the range (s/, s; + dsi) and for the phase factors ■0/ to be on the unit-circle 
segment between and ^; -I- d^i for all I € Bn- In terms of this PDF, taking the averages will involve integration over all 
the real positive sj's and along all the complex unit circles of all i^/'s. 



where the notation /{A^,?/;} means that / depends on all A^'s and all ij^i's in the set {v4^, i/);; ZgBat} (similarly, {s,<^} 
means {si,ipi;l £ Bn}, etc). The full PDF that contains the complete statistical information about the wavefield a(x, i) 
in the infinite x-space can be understood as a large-box limit 



i.e. it is a functional acting on the continuous functions of the wavenumber, Sk and ^k- In the the large box limit there 
is a path-integral version of 



^It is easy to extend the analysis to the infinite Fourier space, kmax = oo. In this case, the full joint PDF would still have to be defined 
as a Af — > oo limit of an A'^-particle PDF, but this limit would have to be taken in such a way that both kmax and the density of the Fourier 
modes tend to infinity simultaneously. 
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The full PDF defined above involves all N modes (for either finite N or in the N ^ oo limit). By integrating out all the 

arguments except for chosen few, one can have reduced statistical distributions. For example, by integrating over all the 
angles and over all but M amplitudes,we have an "M-particle" amplitude PDF, 

which depends only on the M amplitudes marked by labels ji,j2, ■ ■ ■ jJm&Bn. 

Statistical derivations are greatly facilitated by the introduction of a generating functional 

^^^HA,M} = -ip(n e^'^?Vr>> (4) 

where {X,n} = {A(,/x;;Z € Bn} is a set of parameters, XiGTZ and hiGZ. 



^ ' {/*} leBN {n} \ leBN / 



where {/z} = {m G Z]l G Bn} is a set of indices enumerating the angular harmonics and C^^ stands for the inverse 
Laplace transform with respect to all A;. 

2.1 Definition of an essentially RPA field 

A pure RPA fields can be defined as one in which all the phases and amplitudes of the Fourier modes make a set of 
2A'' statistically independent variables and in which all phase factors ijj are uniformly distributed on their respective 
unit circles. In such pure form RPA never survives except for in the un- interesting state of complete thermodynamic 
equilibrium. However, WT closure only requires an approximate RPA which holds up to certain order in small e and 
1/A'' and only in a coarse-grained sense, i.e. for the reduced M-particle objects with M <^ N. Below we give a relaxed 
definition of an (essentially) RPA property which, on one hand, is sufficient for the WT closure and, on the other hand, 
is preserved over the nonlinear time. 

Definition: We will say that the field a is of an essentially RPA type if: 

1. The phase factors are statistically independent and uniformly distributed variables up to O(e^) corrections, i.e. 
where 

7'W«)M = n I'^^'i ) ^^""H^-o, (7) 

is the A^-particle amplitude PDF. In terms of the generating functional 

Z(^){A, n} = Z(^'''){A} n '^(mO [1 + 0{e% (8) 

IGBn 

where 

ZW«){A} = ( n e^'^') = zW{A,m}U=o (9) 
leBN 

is an AT-particle generating function for the amplitude statistics. 
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2. The amplitude variables are independent in a coarse-grained sense, i.e. for each M <^ N modes the M-particle 
amplitude PDF is equal to the product of the one-particle PDF's up to 0{M/N) and o(e^) corrections, 

'Pr^3t.,,M = PnPn ■ ■ ■ PjM [1 + 0{M/N) + O^e')]. (10) 

As a first step in validating the RPA property we will have to prove that the generating functional remains of the form 
over the nonlinear time provided it has this form at t = 0. 

3 Weak-nonlinearity expansion. 

Consider weakly nonlinear dispersive waves in a periodic box with a dispersion relation ujk which allow three-wave 
interactions. Example of such systems include surface capillary waves [21Eli Rossby waves jOj and internal waves in the 
ocean 0. In Fourier space, we have the following Hamiltonian equations, 

oo 
m,n— 1 

where ai = a{ki) is the complex wave amplitude in the interaction representation, ki = 2nl/L is the wavevector, L is 
the box side length, = uJki — ~ '^fcmi '^i = '^fe; is the wave frequency, e ^ 1 is a formal nonlinearity parameter. 
Here, the interaction coefficient VJ^^ is obviously symmetric with respect to m and n but we do not assume any further 
symmetries.^ 

In order to filter out fast oscillations at the wave period, let us seek for the solution at time T such that 27r/tj <C T <C 
1/cije^. The second condition ensures that T is a lot less than the nonlinear evolution time. Now let us use a perturbation 
expansion in small e, 

/rr\ (0) I (1) , 2 (2) /ir,N 

ai[l ) = fl; + eaj + e al . (12) 

Substituting this expansion in (|ll|l we get in the zeroth order af'\T) = ai{0), i.e. the zeroth order term is time 
independent. This corresponds to the fact that in the interaction representation, wave amplitudes are constant in the 
linear approximation. For simplicity, we will write af'\o) = ai, understanding that a quantity is taken at T = if its 
time argument is not mentioned explicitly. The first order is given by 

oc 

a'i'\T) J2 (KLn«n^a„A^„4+„ + 2Vr:a^anArjr+n) , (13) 

m,n— 1 

where Aj^^^ = e^'^^^^^dt — (e*"^"^"^ — l)/iujl^^^^. Iterating one more time we get 

oo 



(14) 



where we introduced E{x, y) — A(x — y)e^y*dt. 

^Some additional symmetries involving permutations of the upper and lower indices arise, e.g., in solids due to the fact that nonlinearity 
is purely due to the potential energy which is a function of the displacement but not the rate of the displacement. Refs. I15lll6lll3l imposed 
such symmetries which immediately rule out the capillary, internal and other waves in fluids for which such properties do not hold. Additional 
symmetries also arise if the action variable is a Fourier transform of a real quantity, e.g., in the Rossby waves 
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4 Evolution of the Generating Functional and Multi-particle PDF 



Let us first derive an evolution equation for the generating functional Z{A,/x} exploiting the separation of the linear and 
nonlinear time scales. ^ To do this, we have to calculate Z at the intermediate time t = T based on its value at i = 0. 
The derivation, although standard for WT, is quite lengthy and will have to be published in a longer paper. Here, we 
will only outline the main steps and give the result. First, we need to substitute the e-expansion of a from (|12|l into the 
expressions e'^^'"^' and ipj^ = ^(-^'^ 1^)'^^ • Second, the phase averaging should be done. Note that, because, we assume 
that initial phase factors are independent at t — with required accuracy, we can do such phase averaging independently 
of the amplitude averaging (which we do not do yet). Thirdly, we take N ^ oo limit followed by T ~ 1/e — > cx) (this 
order of the limits is essential!). Taking into account that lim E{0,x) — T {Trd{x) + iP{-)), and lim |A(x)p = 27rT5(x) 

and, replacing {Z{T) — Z{0))/T by Z (because the nonlinear time ~ T) we have 



S 



Ane' I {(A, + A^^^) |KLP<5«„)5^+„ + 2\Vpj'S{L.J^JS[ 



+2A,- 



s^z 

n 

sz 



1 ()^ 7 

''^\yLS^UJ{^Ln) + I^^P-^^+^-J^™) )dk,dk^dK. (15) 



(5Aj(5A„(5Am 



Here variational derivatives appeared instead of partial derivatives because of the N oo limit. This expression is valid 
up to the [1 + 0(e^)] factor. Equation (|15|l does not contain fj, dependence which means that that these variables separate 
from A's and the solution is a purely-amplitude Z times an arbitrary function of /i's which is going to be stationary in time. 
The latter corresponds to preservation of the initial H(5(/i/) dependence by equation (|15|l which means that no angular 
harmonics of the PDF higher than zeroth will be excited. In the other words, all the phases will remain statistically 
independent and uniformly distributed on with the accuracy of the equation (|15l) integrated over the nonlinear time 
l/e^, i.e. with the 0{e'^) accuracy. This proves the first of the "essential RPA" properties. In fact, this result was already 
obtained before in [IS] for a narrower class of 3- wave systems (see footnote 2). Note that we still have not used any 
assumption about the statistics of A's and, therefore, (|15|l could be used in future for studying systems with random 
phases but correlated amplitudes. 

Taking the inverse Laplace transform of H15|) we have the following equation for the PDF, 

V^-J'^dk„ (16) 

where Fj is a flux of probability in the space of the amplitude Sj , 

= /{(I^^J''5(..^J<5^,+„ + 2|^"„j2<5(c.;.'„)5;+„J.s„s,„g 

|2 r/, .n NX" lT/7 |2 ; 



sr 



2{\V;JH{lu^;JS^^„^ - 2|y^„|2^«J<5i+Js„s„ — } dk^dk^ (17) 



This equation is identical to the Zaslavski-Sagdeev (ZS) ^31 equation (Brout-Prigogine in the physics of crystals context 
|15[ I16|). Note that ZS equation was originally derived in ^3] for a much narrower class of systems, see footnote 2, 
whereas the result above indicates that it is also valid in the most general case of 3-wave systems. Here we should again 
emphasize the importance of the order of limits, — > oo first and e ^ second. Physically this means that the frequency 

^Hereafter we omit superscript (A') in the Af-particle objects if it does not lead to a confusion. 
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resonance is broad enough to cover a great many modes. Some authors, e.g. ZS and BP leave the sum notation in the 
PDF equation even after the e — > hmit taken giving (5(tj"„). One has to be careful interpreting such a formula because 
formally the RHS is null in most of the cases because there may be no exact resonances between the discrete k modes 
(as it is the case, e.g. for the capillary waves). Thus, our functional integral notation is a more accurate way to write the 
result. 



5 In what sense are the amphtudes independent? 

Obviously, the variables sj do not separate in the above equation for the PDF. Substituting 

'p(N,a) _ p(a) p(a) p(a) 

(compare with (|10|) ') into the discrete version of (|17|) we see that it turns into zero on the thermodynamic solution with 
pj"^ = Ljj exp(— WjSj). However, it is not zero for the one-mode PDF pj°^ corresponding to the cascade-type Kolmogorov- 

Zakharov (KZ) spectrum n^^ , i.e. pj"'' = (l/n*^^) exp(— Sj/n*^^) (see next section), nor it is likely to be zero for any other 
PDF of form (|18|l . This means that, even if initially independent, the amplitudes will correlate with each other at the 
nonlinear time. Does this mean that the existing WT theory, and in particular the kinetic equation, is invalid? 

To answer to this question let us differentiate the discrete version of the equation (|15|l with respect to A's to get 
equations for the amplitude moments. We can easily see that 

dt{{AlAl)~{Al){Al))^Oie') (ji,j2 G M (19) 

if {Aj^Aj^AjJ = {A'j^){Aj_^){A'j_^) (with the same accuracy) at t = 0. Similarly, in terms of PDF's 

dt{ptSi'^^''n)-Pn\^n)Pj:\^n))-0{e^) (ji,j2 e S^) (20) 

-^jl,h,j3,j4. ' ; ' *J4) — Pji ' i^jl)Pj2 ' (•^J2)^J3 (^J3 )^J4 (^^4 ) at f = 0. HcrC ^'/j j 3 J 3 J 4 ' ' *J3 ' *i4 ) J Pjl,j2 ' ^h) 

and Pj'^\sj) are the four-particle, two-particle and one-particle PDF's obtained from V by integrating out all but 4,2 or 
1 arguments respectively. One can see that, with accuracy e^, the Fourier modes will remain independent of each other 
in any pair over the nonlinear time if they were independent in every triplet at t = 0. 

Similarly, one can show that the modes will remain independent over the nonlinear time in any subset of M < N 
modes with accuracy Ai/N (and e^) if they were initially independent in every subset of size M + 1. Namely 

PJ^!^:^...,M , : ^.M ) - P^f i^n (^^2 ) • ■ ■ (^,M ) = 0{M/N) + 0{e') 

{ji,j2,---jM e Bn) (21) 

:f p(M+l,a) _ p(a) p(a) p(a) „j- ^ _ n 

The mismatch 0{M/N) arises from some terms in the ZS equation with coinciding indices j. For M = 2 there is only 
one such term in the A^-sum and, therefore, the corresponding error is 0{l/N) which is much less than O(e^) (due to the 
order of the limits in N and e). However, the number of such terms grows as M and the error accumulates to 0{M/N) 
which can greatly exceed O(e^) for sufficiently large M. 

We see that the accuracy with which the modes remain independent in a subset is worse for larger subsets and that 
the independence property is completely lost for subsets approaching in size the entire set, M ^ N . One should not 
worry too much about this loss because N is the biggest parameter in the problem (size of the box) and the modes 
will be independent in all ilZ-subsets no matter how large. Thus, the statistical objects involving any finite number of 
particles are factorisable as products of the one-particle objects and, therefore, the WT theory reduces to considering the 
one-particle objects. This results explains why we re-defined RPA in its relaxed "essential RPA" form. Indeed, in this 
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form RPA is sufficient for the WT cfosure and, on the other hand, it remains vahd over the nonhnear time. In particular, 
only property (|19|) is needed, as far as the amplitude statistics is concerned, for deriving the 3- wave kinetic equation, and 
this fact validates this equation and all of its solutions, including the KZ spectrum which plays an important role in WT. 

The situation where modes can be considered as independent when taken in relatively small sets but should be treated 
as dependent in the context of much larger sets is not so unusual in physics. Consider for example a distribution of 
electrons and ions in plasma. The full A'^-particle distribution function in this case satisfies the Louville equation which is, 
in general, not a separable equation. In other words, the TV-particle distribution function cannot be written as a product 
of N one-particle distribution functions. However, an M-particle distribution can indeed be represented as a product of 
M one-particle distributions if M <C Nu where No is the number of particles in the Debye sphere. We see an interesting 
transition from a an individual to collective behavior when the number of particles approaches Nu. In the special case of 
the one-particle function we have here the famous mean- field Vlasov equation which is valid up to 0{1/Njj) corrections 
(representing particle collisions). 



6 One-particle statistics 

We have established above that the one-point statistics are at the heart of WT theory. All one-point statistical objects 
can be derived from the one-point amplitude generating function. 



which can be obtained from the A''-point Z by taking all /i's and all A's, except for Xj, equal to zero. Substituting such 
values into H15|l we get the following equation for Za, 

^ = A,r;,Z„ + (A,^ry,-A,7.)^, (22) 

where, 

V, = 4^6^ J {\Viij'SlJ{cuiJ + 2\Vjl^f6]lS{iu]i)) nium dhdk^, (23) 

7, = 8ne^ j (\VL?5iJ{oji;)n„, + \V;^\H^Pi5{ujJ){ni ~ n,„)) dhdk^. (24) 



Correspondingly, for the one particle PDF Pa{sj) we have 

dPa dF 



dt ds 



= 0, (25) 





with is a probability flux in the s-space, 

()P 

F = -s,{jPa + v,j^). (26) 

Equations l|22|l and (|25|l where previously obtained and studied in [21] in for four-wave systems. The only difference for 
the four- wave case was different expressions for rj and 7. For the three- wave case, the equation for the PDF was not 
considered before, but equations for its moments were derived and solved in jl9j . In particular, the equation for the first 
moment is nothing but the familiar kinetic equation n = —^n + rj which gives rj = jn for any steady state. This, in turn 
means that in the steady state with _F = we have Pj""^ — (l/rij) exp(— Sj/rij) where rij can be any steady state solution 
of th kinetic equation including the KZ spectrum which plays the central role in WT pi I17j . However, it was shown in 
[21] that there also exist solutions with F ^ which describe WT intermittency. 
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7 Discussion 



In the present paper, we considered the evolution of the full N-particle objects such as the generating functional and the 
probability density function for all the wave amplitudes and their phase factors. We proved that the phase factors, being 
statistically independent and uniform on initially, remain so over the nonlinear evolution time. This result does not 
rely on any assumptions about the amplitude statistics and, therefore, can be used in future for studying systems with 
correlated amplitudes (but random phases). If in addition the initial amplitudes are independent too, then they remain 
so over the nonlinear time in a coarse-grained sense. Namely, all joint PDF's for the number of modes M ^ N split 
into products of the one-particle densities with 0{M /N) accuracy. Thus, the full A^-particle PDF does not get factorized 
as a product of N onc-particle densities and the Fourier modes in the set considered as a whole are not independent. 
However, the wave turbulence closure only deals with the joint objects of the finite size M of variables while taking 
the N ^ oo limit. These objects do get factorized into products and, for the WT purposes, the Fourier modes can be 
interpreted as statistically independent. These results reduce the WT problem to the study of the one-particle amplitude 
PDF's and they validate the generalized RPA technique introduced in Such a study of the one-particle PDF and 

the high-order momenta of the wave amplitudes was done in |^ and the reader is referred to these papers for the 
discussion of WT intermittency. 

Finally, we would like to mention the role of quasi-resonant interactions which, as we saw, do not produce any long- 
term effect at the order considered in this paper. However, these interactions do modify statistics at order as was 
shown in |22| . The correction can be important for the real space correlators which have Gaussian values at the 
order for any (not necessarily Rayleigh) amplitude distributions. 
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